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4^ Abstract 
^— > 

Qh We consider the conformal field theory of N complex massless scalars in 

2 + 1 dimensions, coupled to a U (N) Chern-Simons theory at level k. This 
theory has a 't Hooft large N limit, keeping fixed A = N/k. We compute 
some correlation functions in this theory exactly as a function of A, in the 
OH large N (planar) limit. We show that the results match with the general 

^ predictions of Maldacena and Zhiboedov for the correlators of theories that 

have high-spin symmetries in the large N limit. It has been suggested in the 
past that this theory is dual (in the large N limit) to the Legendre transform 
CN of the theory of fermions coupled to a Chern-Simons gauge field, and our 

results allow us to find the precise mapping between the two theories. We 
find that in the large N limit the theory of N scalars coupled to a U(N)k 
Chern-Simons theory is equivalent to the Legendre transform of the theory 
of k fermions coupled to a U (k) n Chern-Simons theory, thus providing a 
bosonization of the latter theory. We conjecture that perhaps this duality 
is valid also for finite values of N and k, where on the fermionic side we 
should now have (for Nj flavors) a U(k) N _ Nf / 2 theory. Similar results hold 
for real scalars (fermions) coupled to the 0(N)k Chern-Simons theory. 
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1 Introduction and Summary of Results 

Two of the simplest conformal field theories in 2 + 1 dimensions are the theory 
of iV real massless free scalars (p l , in the sector of operators that are singlets of 
O(N), and the "critical O(N) model". The latter may be viewed as the IR limit 
of the deformation of the free theory by g^^ip 1 ) 2 (when fine-tuning the IR scalar 
mass to zero), or equivalently by adding an auxiliary field a and deforming the 
free theory by <j(p l (p l (and fine-tuning a linear term in a to get to a non-trivial 
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fixed point). Due to the latter description we will call the "critical O(N) model" 
the Legendre transform of the free theory 1 . 

Two other simple conformal field theories in 2 + 1 dimensions are the theory 
of N massless free fermions ip a (again limited to singlets of O(N)), and the Gross- 
Neveu model. One can formally define the latter by deforming the free fermion 
theory by (ip a ip a ) 2 , but since this operator is not renormalizable this definition 
does not make much sense (except in the large N limit). As above, one can also 
formally define it by performing a Legendre transform with respect to the operator 
(ip a 'ifj a ) . One way to properly define the Gross-Neveu theory is as a theory that 
has a scalar operator a (analogous to the auxiliary field a above), such that when 
we deform the theory by (<7 2 ) it flows to the free fermion theory (when fine-tuning 
the IR fermion mass to zero). 

The free theories described above both have high-spin symmetries of every 
even spin s = 2,4, 6, • • • , generated by conserved currents J^. In the interacting 
theories these symmetries (for s > 2) are broken, though the effect of the breaking 
is small in the large N limit 2 . In this limit all local operators in these theories are 
products of the high-spin currents and of one additional scalar operator 
(given by (<£>V l ) i n the bosonic case, and by (ijj a ijj a ) in the fermionic case). 

The two bosonic theories are simply related by a Legendre transform with 
respect to J^°\ as are the two fermionic theories, but at first sight there is no 
relation between the bosonic and fermionic theories. As far as we know, the first 
hint for such a relation came by considering their gravity duals. All the theories 
discussed above have a good 1/N expansion (see [1] for a review), so it is natural 
to suggest that they could have classical gravitational duals (by the AdS/CFT 
correspondence [2, 3, -1]) at large N, living on AdS±. These duals should have 
massless high spin fields, to match with the field theory spectrum. Indeed, it was 
suggested in [5, 6] (see also [7, 8, 9]) that the bosonic theories are dual to the type 
A high-spin gravity theories of Vasiliev [10], and the fermionic theories to the type 
B high-spin gravity theories of Vasiliev. In the gravity language, the difference 

1 The two theories are not precisely Legendre transforms, except in the large N limit, since 
we need to flow to the IR and to fine-tune to get to the fixed point. However, we will use this 
name for this type of relation between theories, for lack of a better name. 

2 The divergences of the high-spin currents and their anomalous dimensions are proportional 
to 1/N, but some correlators feel the breaking even in the large N limit. 
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between the two bosonic (fermionic) theories is just a different choice of boundary 
conditions for the bulk scalar field which is dual to J(°>] in the classical gravity 
limit this is equivalent to a Legendre transform [11, 12, 13]. Strong evidence for 
this equivalence was found in [14, 15], and suggested derivations of the equivalence 
were presented in [16, 17, 18, 19]. 

The type A and type B high-spin gravity theories mentioned above have the 
same spectrum but different interactions. However, there is a family of high-spin 
gravity theories labeled by a parameter 9 (appearing in the Vasiliev equations of 
motion) that interpolates between the type A and type B theories, such that the 
bosonic theories arise when 9 = and the fermionic ones when 9 = ir/2. This 
hints that perhaps there is a family of field theories (that is continuous, at least in 
the large N limit) which interpolates between the bosonic and fermionic theories 
described above. The theories with 9 ^ 0,7r/2 are not parity- invariant. It was 
conjectured in [ ] that they arise by coupling the bosonic/fermionic theories to 
O(N) Chern-Simons theories at level k. For infinite k this just implements the 
reduction to singlets of O(N), but for finite A = N/k in the large N limit it 
provides a parity-breaking modification of the bosonic and fermionic theories, by 
a parameter that is continuous in the large N limit (though it is discrete for finite 
N). This conjecture was tested in [20, 21], where it was shown that the spectrum 
of operators in these theories is independent of A in the large N limit. This is 
consistent with the fact that the high-spin fields in the dual gravity theory are 
classically massless; presumably they acquire a mass (whenever the field theory is 
not free) at order 1/N by loop corrections in the bulk [22]. 

Significant support for this conjectured relation between the bosonic and fermionic 
theories was recently given in [23, 24], where it was shown that the large N cor- 
relation functions in these theories could be computed just by knowing that they 
have a high-spin symmetry that is broken by 1/N effects (plus some additional 
technical assumptions). Maldacena and Zhiboedov showed that these correlation 
functions could be expressed in terms of two (or three) parameters appearing in 
the non-conservation equations of the high-spin currents. It is natural to assume 
that these parameters map to the parameters N and k mentioned above, but the 
precise relation is not known (when there is a third parameter it can be identified 
with the coefficient of an extra ge((p % (f 1 ) 3 interaction, that is exactly marginal in 
the large N limit [ ]). Maldacena and Zhiboedov found that if one starts from the 



4 



free bosonic theory and turns on the coupling (in their language, turning on the 
breaking of the high-spin symmetry) , then in the limit of infinite coupling one ends 
up with the correlation functions of the Gross-Neveu model (the Legendre trans- 
form of the free fermionic theory). The Legendre-transformed statement is also 
true: by starting from the correlation functions of the free fermionic theory and 
turning on the coupling, they found that in the limit of infinite coupling one ends 
up with the "critical O(N) model" (the Legendre transform of the free bosonic 
theory). 

This suggests that for each theory mentioned above, there are actually four 
different ways to describe it and to compute its large iV correlation functions : 
(a) It can be described as a theory of N s massless scalars coupled to a 0(N s )k a 
Chern-Simons theory; (b) It can be described as a theory of iVf er . massless fermions 
coupled to a 0(N{ ei ,)] Cfei Chern-Simons theory; (c) It can be described purely 
algebraically as a theory with a slightly-broken high-spin symmetry, and its large 
iV correlators can be expressed in terms of the parameters N and A of [24]; (d) It 
can be described as Vasiliev's high-spin gravity theory with a parameter 9. The 
first two descriptions exist also for finite N, while in the latter two it is only known 
how to compute in the planar limit. The discussion above implies that in the large 
N limit all of these descriptions are equivalent (up to possible Legendre transforms, 
in particular between the first and second descriptions), but the precise mapping 
between them is not yet known. 

Our main goal in this paper is to clarify the relation between the first three 
descriptions of these theories, and to compute the precise mapping between their 
parameters (in the large N limit) 3 . All the O(N) theories discussed above have also 
U (N) versions, where we start with iV complex scalars (fermions) and couple them 
to a U(N)k Chern-Simons theory. The U(N) theories have conserved currents 
(in the large N limit) also for odd spins s, but in the large N limit the correlators 
of the even spin operators in these theories are equivalent to those of the O(N) 
theories. In this paper we will perform computations for the U (N) case, since in 
this case there is a spin-one current that takes a simpler form than the higher spin 
currents; however our conclusions should be equally valid for the O(N) case. 

3 We will not discuss here the mapping to the gravitational side. The results of [20] for even 
contributions to correlators imply that in this mapping A = tan(0), and this is confirmed by 
more detailed computations in [25]. 
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Most of our paper is devoted to performing exact computations of large N 
correlation functions in the bosonic theories described above. After introducing 
our theories and methods in section 2, we compute in section 3 some exact planar 
correlation functions of the scalar fields in these theories. These are then used 
in section 4 to compute exactly (as a function of A) some planar correlators of 
gauge-invariant operators in the bosonic theories. Specifically we compute some 
2-point and 3-point functions of the operators with s — 0, 1, 2. 

Section 5 contains our main results. We begin by matching our results for the 
theories of scalars coupled to a U(N)k Chern-Simons theory to the general results 
of [24]. There are two different standard definitions of the Chern-Simons level. We 
will denote by k the definition we will use in our computations (which is sometimes 
called the "renormalized coupling"; it arises, for instance, from regularization by 
dimensional reduction). However, defining the level by a regularization using a 
Yang-Mills term at high energies gives a different (at one-loop) definition of the 
level, &ym = k — N. Some of our results are simpler to express in the first lan- 
guage and some in the second language; of course the translation between them is 
straightforward, and we will try to carefully distinguish the two everywhere. The 
matching to [24] works better using k and A = N/k; in these variables we find that 
the parameters of [24] are related to the bosonic rank and level by 4 



In particular A diverges as A — > 1, which corresponds to an infinite coupling 
(N/kyyi — > oo) using the Yang-Mills regularization, and in this limit we should 
approach the fermionic theory. 

The mapping to the fermionic theory is nicer to describe using the Yang-Mills 
regularization, which we will use in the rest of this introduction 5 . Taking the strong 
coupling limit, we find that the theory of iV scalars coupled to a U(N) kYM Chern- 
Simons theory matches (in the large N,kyM limit) to the (Legendre transform of 
the) theory of kyM fermions coupled to a U(kyM)N Chern-Simons theory 6 . The 

4 The same value of A is found in supersymmetric generalizations of these theories, by com- 
pletely different methods, in [25]. 

5 We will always use this convention whenever we have kyM appearing, either in the level or 
in the rank of our theories. 

6 In our computations we do not determine the sign of the Chern-Simons coupling, so the level 
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two Chern-Simons theories in question are related by level-rank duality [26, 27, 28], 
and our claim is that coupling one of them to a massless scalar in the fundamental 
representation is exactly equivalent (up to a Legendre transform) to coupling the 
other to a massless fermion in the fundamental representation. This may be viewed 
as a bosonization of the fermionic theory, expressing it purely in bosonic variables. 
The generalization to the O(N) case is straightforward, mapping the theory of 
N real scalars coupled to an 0(N)k YM Chern-Simons theory to the (Legendre 
transform of the) theory of &ym real fermions coupled to an 0(/cym)at Chern- 
Simons theory. 

Our results show that the bosonic and fermionic theories have the same corre- 
lation functions at large N, and thus provide strong evidence for their equivalence 
in this limit. More precisely, we derive this particular matching of parameters just 
in the strong coupling limit of the bosonic theory, which is the weak coupling limit 
of the fermionic theory, but it seems natural that it should extend throughout the 
parameter space; this is confirmed by preliminary computations of exact 2-point 
and 3-point correlators in the fermionic theory [29] . The mapping of the previous 
paragraph contradicts computations of the thermal free energy presented in [20], 
which is why it was not found already in [20, 24]; we discuss in the final section 
why these computations do not give the correct answers for the free energy at 
finite A. 

While so far we only discussed the large N limit, it is natural to conjecture 
that perhaps the equivalence between the two theories is valid also at finite N (this 
is true for the level-rank duality, and also in similar dualities for supersymmetric 
theories [30, 31]). If this is not correct, then the scalar and fermion theories 
would provide two different quantum generalizations of the same classical high-spin 
gravity theory. Unfortunately, it is hard to test this conjecture, since at least one 
side of the duality is always strongly coupled, and we do not know how to perform 
exact computations at finite N. All we have so far is a weak test of this finite N 
duality, by comparing the mass deformations on both sides. This test suggests that 
at finite N the precise mapping is from a U (N)k YM bosonic theory to a U (&ym) jv-i/2 
fermionic theory. If we have Nf flavors in the fundamental representation, a natural 
generalization would map the theory of Nf massless scalars in the fundamental of 
U (N) coupled to the U {N)k YM Chern-Simons theory to the (Legendre transform of 

may also be (-N). 
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the) theory of Nf massless fermions in the fundamental of £/(/cym)j coupled to the 
U (&ym) N—Nf/2 Chern-Simons theory. It would be interesting to find ways to test 
this conjecture. In particular it may be interesting to see how the scalars become 
fermions and vice versa; one can try to analyze this 7 by looking at open Wilson 
lines ending on scalars/fermions and checking which anyonic statistics the ends of 
these Wilson lines obey, function of N and k YM . 

2 Vector Model with Chern-Simons Interactions 

Consider the theory of a complex scalar field <p in the fundamental representation 
of U(N), coupled to gauge bosons with Chern-Simons interactions at level k 
in three Euclidean dimensions. The action is 

S = Scs + S 

scalar ■ (2) 

The Chern-Simons action Sqs is given by 

S cs = ^ / Tr N ^AAdA+ 2 -A^ = J d*x (a^A; + \f^A-J^j . 

(3) 

The trace is taken in the fundamental representation, and we normalize the gener- 
ators with C\(N) = — 1/2 8 . With this normalization the theory is gauge invariant 
if k G Z [32, 33]. The scalar field action is 

Scalar = J d 3 X (\V ^ + ^(</>V) 3 ) > ( 4 ) 

where = <9 M + A^. 

We work in the 't Hooft large N limit, keeping A = ^ and A6 fixed. In this 
limit, the theory (2) is conformal to all orders in perturbation theory in A and A6 
21]. In the planar limit, our theory is closely related to the 0(N) theory of N 
real scalar fields coupled to an 0(N) Chern-Simons theory at level k; in this limit 
the latter theory (with N = 2N) is simply a projection of our theory, keeping only 

7 We thank N. Itzhaki for this suggestion. 

8 For a representation R of U(N) with generators {T a }, C\(R) is defined by 
Tr R (T a T b ) = d(R)S ab . 
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some of its operators. For finite N the two theories are not equivalent, but all 
of our large N computations (except for the ones that involve operators that are 
projected out when going to O(N)) can easily be translated into computations in 
the O(N) theory as well. 

Let us define light-cone coordinates by 1 ± ix 2 ) I \/2. We work in 

light-cone gauge, A__ = 9 . With this choice of gauge, the A AAA A self-interaction 
of the gauge field vanishes, and the seagull term (f^A 2 ^ is also simplified as we shall 
see. This greatly reduces the complexity of perturbative calculations. The utility 
of light-cone gauge in theories of this kind was first noticed in [20] . 

In this gauge the gluon propagator is 

(A;(-p)Al(q)) = G u ,{p)5 ab ■ (27r) 3 5 3 (g - p) , 

Ani 1 
k p + 

and the other components of G Vil vanish. At leading order in the 't Hooft large N 
limit, where scalar loops can be ignored, the gluon propagator in this gauge does 
not receive corrections, because the A A A A A interaction vanishes and the ghosts 
are decoupled. 

To regulate the theory, we use dimensional regularization in the direction x 3 , 
and a cutoff A on the momentum in the 1-2 plane. The cutoff regulator breaks 
Lorentz invariance (as does our choice of gauge), conformal invariance and gauge 
invariance. However, we only encounter power-law divergences, and the counter- 
terms used to subtract those divergences are completely fixed by demanding that 
the continuum theory is conformally invariant (in general these counter-terms are 
not gauge-invariant or Lorentz-invariant). The fact that our final results for the 
correlation functions of gauge-invariant operators are consistent with the analysis 
of [24] (see section 5) gives strong evidence that Lorentz and gauge invariance are 
restored in the continuum limit. 

The value of k in Chern-Simons theories depends on the regularization; as in 
any other theory, coupling constants in different regularizations are not the same 
at higher orders in perturbation theory. In particular, as shown in [34], the value 
of k in our regularization (which is sometimes called the "renormalized coupling" ) 

9 In Euclidean space A + — A_, but we keep A + ^ 0. One can think of this prescription as an 
analytic continuation of light-cone gauge in Minkowski space. 
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differs from the value of k using a regularization involving a Yang-Mills term in 
the UV by A; us = k Y M + N = k YM (l + Aym) (where we define Aym = N/k YM ) w . 
In terms of the 't Hooft coupling we have A us = Aym/(1 + Aym), so that the two 
couplings agree perturbatively, but the maximal value of A that can be achieved 
in the Yang-Mills regularization is A = 1 (which corresponds to Aym — ?■ oo, or 
k Y M -> 0). 

In the large N limit, the spectrum of operators of the theory (2) includes a 
single primary operator for each integer spin s > 0, with conformal dimension 
A = s+ l + 0(l/N) [ ]. Each can be written as a symmetric, traceless tensor 
that is schematically given by 

41m s =^---^> 4 + ---- ( 6 ) 

All other primaries are products of these "single-trace" operators. The currents 
J M = jjp , T^y = jffl correspond to the unbroken U(l) and Poincare symmetries 
and are conserved 11 . The currents with s > 2 are generally not conserved when 
A 7^ 0. In this work we will need the explicit form of the following operators : 



T 



+ <W('")- (7) 



Here 2^ = d ^ l + A fJi , ^ = <9 M — A^, and parentheses around indices denote aver- 
aging over symmetric permutations. The trace terms in T^ u will not be important 
for us in this paper. The currents in (7) are canonically normalized, namely, the 
charges Q = J d 2 xJ , = J d 2 xT 0f , obey [Q, \ = <p\ [P M , 0*] = -id^\ 

For completeness we now review the results of [24] which will be needed in 
section 5. 



10 Without loss of generality we assume k > 0. A parity transformation takes k — > —k, A — > —A. 

11 Note that the naive U(l) global symmetry acting on the complex scalar fields is gauged, 
but the equations of motion of the gauge field imply that the symmetry generated by the 
topologically-conserved current — *tr((L4), acts on the scalar fields in the same way as 
the naive U(l) global symmetry, up to an overall factor. 
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2.1 Review of the Results of Maldacena and Zhiboedov 

Our theory belongs, in the large N limit, to the general class of conformal field 
theories studied in [24], which were named "quasi-boson" theories. In general, 
a CFT belongs to this class of theories if it has a large N expansion, and its 
large N spectrum of operators includes conserved high-spin currents of even spins 
and a dimension one scalar operator. The large N expansion parameter, denoted 
by N, was defined in [2 ] to be proportional to the two-point function of the 
energy-momentum tensor, and the proportionality constant was fixed by requiring 
that N — 1 for a free real boson. The theory described in this section, and the 
Gross-Neveu model coupled to U(N)k Chern-Simons interactions, are particular 
examples of "quasi-boson" theories. 

A closely related class of theories, called "quasi-fermion" theories, was defined 
in [24] by the same properties that were used to define the "quasi-boson" theories, 
except that the scalar primary has dimension 2 + 0(1/N). Following [24] we 
will denote this scalar by to avoid confusion with the quasi-boson case. The 
Legendre transform of the theory described in this section, and the theory of N 
free fermions coupled to U(N)k Chern-Simons interactions, are both examples of 
"quasi-fermion" theories. 

The conservation of high-spin currents in either the quasi-boson or quasi- 
fermion theories is generally violated by double-trace and triple-trace operators. 
For example, the non-conservation of the s = 4 current in the quasi-boson theory 
takes the schematic form (with appropriate derivatives on the right-hand side) 

d ■ JW = a 2 J^J^ + a 3 ( J (0) J {0) J (0) + J (2) J (0) J (0) ) , (8) 

and a similar equation containing only the double-trace term applies in the quasi- 
fermion case. Using (8) Maldacena and Zhiboedov derived a sequence of Ward 
identities, which enabled them to express all the planar 3-point functions in terms 
of N, A oc Na 2 and also in terms of a% in the quasi-boson case. 

The 3-point functions of single-trace primaries in the quasi-boson and quasi- 
fermion theories are constrained to contain at most three different conformal struc- 
tures in the large N limit [24], 

+ j slS2S3 (J {si) J {s2) J {s3) Ud- (9) 
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In the above equation (-)bos. and (-)fer. refer to the correlators in the theory of a real 
free boson and fermion respectively, while (-)odd refers to independent conformally- 
invariant contributions to the correlators that appear only in interacting theories. 
Explicit expressions for these structures can be found in [35, 36, 23]. 

To extract meaningful information from the 3-point functions (9) one has to 
specify the precise normalization of the operators. In [ ! I] a normalization was 
chosen such that the 2-point functions of currents are equal in the quasi-boson 
and quasi- fermion theories, 

(J«(x)J< a >(0)> = N(J (s) (x)J^(0)) hos . , 0. (10) 

The 2-point functions of the scalar operators were chosen to be 

(jV>{z) J< \0)) = -^(J(°)(x) J(°)(0)) bo , , (11) 
1 + \ z 

(jV>{x) J (0) (0)> = -X{jW>(x) J (o) (0))fe, . (12) 

With these normalizations the coefficients of the conformal structures in the 
3-point functions (9) in the quasi-boson case were found to be (in the large N 
limit) 

1 A 2 A 

1 + A z 1 + \ z 1 + \ z 

a SlS2 o = N 1 , 7 SlS2 o = N X , (13) 

1 + A z 1 + \ z 

a Sl oo = N 



1 + A 2 

* 1 f N V 

where z is a specific constant (denoted by zi{rtQ ec boson ) 3 in [24]). 
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Similarly, for the quasi-fermion theories 

— n r — at \ — m A q f 

1 + A qf 1 + A qf 1 + A qf 



^o = ^, l slS2 o = N-^ } (14) 
1 + A qf 1 + A qf 



/3 Sl oo = ^ ' 



1 + 
7ooo = • 

The coefficients which are not specified in the above equations correspond to struc- 
tures which are inconsistent with conformal invariance. 

To obtain the correlators of the critical bosonic vector model with Chern- 
Simons interactions it is convenient to define J^ t bos = A q fJ(°). In terms of 
Jcrit bos the 2-point function (12) and /3 Si qo are multiplied by XL , while /3 SlS2 g and 
7 SiS2 5 are multiplied by A q f. The correlators of the critical model with A = are 
then obtained by taking the A q f — > oo limit in (12) and (14). The correlators of 
the critical fermionic vector model are obtained from (11) and (13) in the same 
way. 

The normalization of we use in this paper (7) is different than the one 
in (11), and the precise relation between the two normalizations will be fixed in 
section 5. Note that in [ ] only correlation functions of even-spin operators (that 
appear both in the O(N) and in the U(N) theory) were computed, while many of 
our results below involve the spin-one current However, since the odd-spin 

currents are part of the same high-spin symmetry algebra as the even-spin currents, 
and the results of [24] (10), (13), (14) are (as far as the coupling-dependence goes) 
independent of the spins for all s > 0, we assume that the same results hold also 
for odd-spin correlation functions. Our results will provide a consistency check on 
this natural assumption. 

3 Exact Scalar Correlators 

In this section we compute 2-point and 4-point scalar correlators exactly in pla- 
nar perturbation theory. These correlation functions are not gauge-invariant, but 
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they will be useful for computing gauge-invariant correlation functions in the next 
section. The computations are shown in some detail, to illustrate the techniques 
that will be useful later on. 

3.1 Scalar Propagator 

Let us denote the self-energy of the scalar field by S(p;A), such that the full 
propagator is 

(4(pW(i)) = p2 _ A) ■ (2k) 3 5 3 (p+ q) . (15) 

We will compute the self-energy in light-cone gauge at large N by solving its 
bootstrap equation, shown in figure 1. Only non- vanishing planar diagrams are 
shown. The diagram of order A, with a single gluon line, vanishes due to parity 12 . 
Diagrams that involve the A A AAA vertex vanish in light-cone gauge. Diagrams in 
which two seagull vertices are connected by a gluon line also vanish; this is because 
the seagull term 0' A^A^cf) reduces in this gauge to (p^A^A^, and the only non- 
vanishing gluon propagator is (A 3 A + ). All other diagrams which do not appear in 
figure 1, including those involving a (|0| 2 ) 3 vertex, are subleading in 1/N. 




Figure 1: Bootstrap equation for the scalar self-energy. A filled circle denotes the 
full scalar propagator. 

Let us solve the equation of figure 1. The diagrams on the right-hand side of 

12 While the Chern-Simons interactions break parity, our theory is invariant under the combined 
operation of parity plus A t— > —A. Corrections to the self-energy that come with an odd power of 
A must therefore be accompanied by an e^ up symbol, but in S(p; A) there is only a single external 
momentum to saturate it. 
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the bootstrap equation are given by 

d d k d d l {l+pYG^{l-p)G v Ak-p){k + P y 



-(T a {T a ,T b }T b ) ji J 



{2n) d {27r) d [P-S(fc)][/ 2 -S(/)] 

, d d k d d l (I + vYG u Jl - p)G p '(k - l)(k + l) p 
(b) + (c) = 2({T a ,T b }T a T b ) Jl 1 V FJ ^ y) pK 



(27r) rf (27r) d [k 2 - E(Jfc)] [I 2 - £(/)] 



Using the gluon propagator (5), in the planar limit we have 

d d k d d l (I +p) + (k +p) + 1 



(16) 



(a) = -4(ttA)^ 



(2ir) d (27r) d (I - p) + {k - p)+ [k 2 - S(A;)][/ 2 - £( 



. . , . , f d d k d d l (l+p) + (k + l) + 1 . . 

(b) + (c) = H«Xf5> J p^^ l, _ P J )+ \ k _ / )+ [k2 _ m][l2 _ m ■ (17) 

Notice that light-cone gauge preserves the symmetry of rotations in the 1-2 plane, 



which act as A_ e ld A_. Hence, S(p) = E(p s ,p 3 ), where p s = \/p\ + p\. Fur- 
ther, the expressions (17) are independent of ps, and therefore we can guess that 
S(p; A) = f(X)pl 

Let us now focus on diagram (a). The integrals over the p% component of the 
loop momenta in (17) can be carried out, 

d 1 -'^ 1 f d l ~ e k^ 1 1 



(2ny-< k 2 - Z(k) J (27ry^(k 3 ) 2 + (l-f(X))k 2 2k s \/\ - /(A) ' 

(18) 

where we have taken e — > at the end. For the remaining integrals in the 1-2 
momentum plane we use polar coordinates, with k^ = k s e ±ldk / \/2. The angular 
integral can be treated as a contour integral, for example 

27r ,„k + +p + f dz z + V2p + /k s r . .. 

d6 k — P — = & -. J/ n = 2vr[20(A; s - p s ) - 1] , (19) 



o 



k+ - p + J iz z - y/2p+/k s 



where z = e l6k is integrated over the unit circle and is the step function. The 
remaining radial integrals (performed up to the cutoff A) are trivial, and the result 
is 

-(2j9 s -A) 2 , (b) + (c)= - XH } ( ^ p s (p s -k). (20) 



4(1 - /(A)) v " ' ' v " v ' l-f(\y 
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(p,k,q;\) 



n 

P + Q t 



+ q T 

ii 



'2 

A P 

i k 

h 



Figure 2: Connected diagrams in scalar 4-point function. 

The self-energy (a) + (6) + (c) is therefore a pure divergence, 

A 2 

A 2 . 



4(1 -/(A)) 



(21) 



This divergence is subtracted with a mass counterterm (f)'(f), which is determined 
uniquely by conformal invariance. In the continuum theory the result is therefore 
X = 0, namely the scalar propagator in our gauge does not receive any corrections. 



3.2 Scalar 4-Point Function 

In this section we compute the connected scalar 4-point function 

(<f> H (P + Q)4>U-p)Ai- k - <lW 2 ( k ))i ( 22 ) 

shown in figure 2. From now on the overall factor (27r) 3 <5 3 (^ k) is implicit in all 
correlators we write. Without loss of generality we focus on the terms proportional 
to Si^dj^ (the other terms are related to this by permutations of the momenta). 

Let us first consider the sum of 1-loop diagrams that contribute to W (the 
correlator (22) without the delta function) and that include an A 2 (p 2 vertex. There 
are six such diagrams: the ones shown in figure 3, and their reflections along the 
vertical axis. It is easy to compute these diagrams using the methods of the 
previous section; their sum is 

2nNX2 { (k 2 - + p)+ [ik + q)s ~ °° + q)s ~ h + Ps] ~ A } ' (23) 
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Figure 3: 1-loop diagrams that include the A 2 (p 2 vertex, up to reflections. 



p + q 

k + q 

31 



'2 
A P 

i k 

32 



Figure 4: Bootstrap equation for the connected 4-point function, when q = 0. 



The linear divergence is subtracted by a (</>V) 2 counterterm, which is uniquely 
determined by conformal invariance. From now on we will take q ± — for sim- 
plicity, even though this is not the most general 4-point function. For this choice 
of momentum the seagull vertex does not contribute to the 4-point function, and 
W becomes a sum of ladder diagrams. The bootstrap equation for the 4-point 
function with q + = is shown in figure 4. In writing it we used the fact that the 
scalar propagator does not receive quantum corrections in our gauge. 

Let us now restrict the form of Wi^j^ij), k, q; A). Using dimensional analysis 
we can write 



W ili2jlj2 (p,k,q;X) 



~ (k 



Pa 

V<?3 13 



;A 



(24) 



By computing the first few contributions to W explicitly, we can see that they 
depend only on the variables 



x 



k s 


Ps 






' y = n ' 


z 


N 


m\ 





(k + pY 
(k-p)- 



A' 



A 

M 



A = A sign(g 3 



(25) 
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and that the dependence on z is at most linear 13 . Let us assume that this holds 
to all orders, and write 



W{x, y, z; A) = W (x, y,X) + z W^x, y; A) 



(26) 



Returning to the bootstrap equation, after carrying out the loop integrals over the 
3-component and the angular directions, the equation can be written as 

W + W x z = -tm-^z + 2iX(I + h) , 



h = N 



k 3 JO 
A rp 3 

+2z 

Jk s 



dk 



4k? + q! ' 
4k? + qj 



2 > 



(27) 



where x' = r^-. Let us first solve for W\ by equating the coefficients of z and 
differentiating with respect to k s . The result is 

W x {x, y) = C7 1 ( 2/ ) e - 2darctan ( 2 ^ . (28) 

Here Ci(y) is an integration "constant" that is determined, by plugging W\ back 
into (27), to be 



CM 



4mX 



: 2i\ arctan(2y) 



(29) 



One can now solve for Wo similarly, and plug the result in (24). We find 

4vrA|g 3 | 



(\q\,x,y,z; A) 



g ± =0 



"ili2 "3132 



N 



tan ( A arctan (2A'^ 



iz 



exp 



2iA (arctan(2y) — arctan(2a;)) 



(30) 



13 We choose to work with parity-invariant variables. Up to order A 3 , we find 
W\ oW 



. .A 

-47TI Z , 

N ' 

A 2 



W^loCA 2 ) = 47r — [arctan(2A') + 2z (arctan(2y) — arctan(2x))] , 

A 3 r 2 
^^(a 3 ) = ^ 7r *^y arc t an (2A') (arctan(2y) — arctan(2x)) + z (arctan(2y) — arctan(2a;)) 



where x, y, z are defined above. 
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Figure 5: The vertex (J^'cjxft). A cross denotes a J(°> insertion in the free theory, 
and a circled cross denotes the exact vertex. The hatched ellipse denotes diagrams 
in which the 4 scalar lines are connected. 

4 Exact Gauge-Invariant Correlators 

In this section we compute several exact 2-point and 3-point functions of single- 
trace primary operators in momentum space. Using the results of section 3 one 
can compute any such correlator, with all external momenta pointing in the 3- 
direction, by computing a finite number of integrals. A useful step is to first 
compute exact vertices of the form (J^(—q)(j)(f)^), since they encode all the infor- 
mation of the scalar 4-point function inside 2-point and 3-point gauge-invariant 
correlators. We begin by computing correlators of = <p^(p, and then generalize 
to other operators. 



4.1 /J (0) J (0) 



Let us compute the J(°> vertex 

(J®(-q)$(k)P) = 5}V ( q ,k; A), (31) 

for the special case of q ± = 0, where we can use the results of the previous section. 
(From here on we do not always explicitly write the momentum of the last operator 
in our correlation functions, which is fixed by momentum conservation.) As shown 
in figure 5, we can write this vertex as a sum of a free piece, and a piece from the 
connected 4-scalar function W. The free piece is just <5*. The connected piece is 
given by (using the variables of (25)) 



N\qM 



d d p W(x, y, z; A) 
(27r) d p 2 (p + q) 2 



2g— 2iA arctan(2x) 
J _|_ g— 2iAarctan(2A') 



(32) 
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Figure 6: The diagrams contributing to (J^(-q)J^). 



where the integral can be computed by the same method as before. Crucially, the 
radial integration can be carried out analytically. Adding the free piece, we find 

o 2iX arctan(2x) 

V (x; A) = = . (33) 

v ' ' l _|_ e -2iAarctan(2A') v ' 

It is now easy to compute the J^ ' 2-point function, shown in figure 6. Note 
that replacing a single insertion by the exact vertex (33) accounts for all the 
diagrams without any double-counting. Thus, 

<^W>>^/^^^^, ,34) 
J (2iT) d k 2 (k + q) 2 4 |g 3 | n \ 

where in writing the final result we took A — > oo. Note that the result is even 
under A — > —A, as expected. The obvious Lorentz-invariant generalization of this 
result to any q is 

</o W) (35) 

4 \q\ 7rA 

The momentum dependence of (35) is determined by conformal invariance; in 
position space this translates into (j(°'(x)j(°'(y)) oc l/(x — y) 2 . The dependence 
on A is not fixed by the symmetries, and we observe that the 2-point function 
diverges at A = 1, and becomes negative afterwards, suggesting that our theories 
only make sense up to A = 1. As mentioned above, if we define the Chern-Simons 
level using a regularization involving a Yang-Mills theory, the value A = 1 is 
actually the maximal allowed value, so this is not too surprising (similar results 
were found also for fermions coupled to Chern-Simons theory in [20]). 

4.2 (jWjWj® 



Using the exact vertex (33) we can compute the 3-point function of scalar op- 
erators ( J(°)(— g) g') q")}, with all momenta q, q', q" in the x 3 -direction. 
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(a) (b) (c) (d) (e) 

Figure 7: The diagrams contributing to (J^°\— q) q') J^\— q")). 

Using the vertices defined above, all contributions to this 3-point function are in- 
cluded in the diagrams of figure 7. All the diagrams turn out to be finite, so we 
remove the cutoff in the expressions below. Diagram (a) evaluates to 
m , d 3 k V (q,k)V (q',k + q')V (q",k~q") 



(2tt) 3 k 2 (k + q) 2 (k - q") 2 

N tan (f) 1 



, , . (36) 
2ttAcos 2 (f ) \q\\q'\\q"\ 

Summing over diagrams (b), (c) and (d), and integrating over the 3 and angular 
directions gives 

(b) + (c) + (d) = -*» r^kjup. 



vr J ^ (Ak 2 + q 2 )(Al 2 + q' 2 ) (Ap 2 s + q" 2 ) 

{{2Q{k B -p s )-l){2Q{l s -p s )-l) 
+ (29(p s - k s ) - 1) (20(/ s - k s ) - 1) 
+ (2Q(k s - l s ) - 1) (26(p fl - l s ) - 1) } 



vr Vio 4fc2 + q 2 J V7o 4/2 + q' 2 J \J ^ Ap 2 s + q" 2 
iVtan 3 (f) 1 



87rA loi^'l |<?"| 

where in the second equality we used the fact that 

(2e{k s - p s ) - 1) {2Q{l s - p s ) - 1) 
+ (2&( Ps - k s ) - 1) (20(/ s - k s ) - 1) 

+ (20(A; S - I.) - 1) (2e( Ps - l„) - 1) = 1 . (38) 
The remaining diagram (e) contributes if Xq ^ 0. It gives 

(e) = NX ^(f) 1 (39) 
{ ) 64vr 3 A 3 \q\\q'\\q"\ ' 1 J 
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Note that (J^j^jw) is the only 3-point function which receives a contribution 
from the (</>V) 3 vertex in the planar limit. This is since diagram (e), when com- 
puted in momentum space, factorizes into a product of three 2-point functions of 
J(°) with the inserted operators: (J^J^J^) (c) oc (J^ J®)(jW j(°))(j( fl ") 
and these 2-point functions vanish by conformal invariance unless s = s' = s" = 0. 
Summing up all the contributions we obtain 

(jW(- q )jW(- q ')j(°\- q ")) = - 

z 

(40) 

which has the correct momentum dependence required by conformal invariance (it 
is uniquely determined by the result we computed for q ± = q f± = q n± = 0). 



ta"(f) 



tarn 



cos- 



7rA 



in 



2 A 2 



\q\ \q'\ \q" 



4.3 (jWjW) 

Next, we compute two correlators that involve J^>. 

In this subsection we compute («/_(— q)J+), with = 14 . For this purpose 
let us first compute the J_ vertex, namely (J-(— q)<jy<j>) with g = 0. In fact, it is 
no more work to compute the J vertex 

(J_..._(-g)0 J t (A;)0 i ) = ^(g,A;;A), (41) 

for any current with s > 0, as long as q ± = 0. Indeed, in our gauge the 
bootstrap equation for this has the same form as the equation for the vertex 
(figure 5), since A_ = so there are no additional diagrams where a gauge field 
is connected to J^ s \ Further, by Lorentz invariance the free piece in (41) is given 
(when q* 1 = 0) by a s (k + ) s , where the factor a s depends on the normalization of 
the current. In computing the connected piece of the bootstrap equation, the only 
modification to (32) is to insert a s (p + ) s in the integral. The computation carries 
through, and we find for all s > 1 

V = a ^+^ e 2 *^( arctan ( 2A ')- arctan ( 2:r )) ^42) 

14 Note that when q — 0, current conservation implies 93 (^(9) ■ ■ • ) = up to contact terms, 
and therefore correlators involving the J3 component of do not contain information that is 
useful for our purposes. 
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Figure 8: The diagrams contributing to («/-(— q)J+)- 

U^k-X) = 1 — B 3 = X + 

k J + 



— + X d — 13 + 

4 reflections . 

Figure 9: Diagrams of (J + 0^0) in which a gauge field is connected to J + . 

Note that a\ — 2 for the canonically-normalized 

The 2-point function (J-(— q)J+) receives the contributions shown in figure 8, 
where the circled insertion is the exact vertex given by V\ of (42), and the squared 
insertion is defined in figure 9: it accounts for the diagrams in which a gauge field 
A + is connected to J + . The squared insertion of figure 9 evaluates to 

U l{q3 ,t,X)= k - + f^ k - 6i. (43) 

The exact 2-point function (J_(— q)J+) is then given by 

d d p V 1 (q,p)U 1 (-q,p + q) 



(J-{-q)J+) = N 



(2-K) d p 2 (p + qf 



16 7rA 4tt ( ' 

In the second line we took A — > oo, carefully keeping divergent terms. 

This result has two peculiar features. First, the linear divergence is a contact 
term that violates conformal invariance. If we denote the background gauge field 
that couples to J M by ^4 M , then this divergence can be subtracted by a mass coun- 
terterm A^A^. The second peculiar feature of (44) is the appearance of a A-odd 
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part, equal to 

iNq 3 cos(ttA) - 1 

~6 7rA • (45) 

This piece violates parity, while {J J) is parity-even in a conformal theory [37, 
35]. However, this piece can come from a contact term {J^(q)J u ) ~ e Ml/p g p that 
is conformally-invariant, and gives a derivative of a delta function in position 
space. It corresponds to the appearance of a Chern-Simons term ^ J A A dA m 
the generating functional F[A, ...]. In many cases, contact terms are scheme- 
dependent and therefore do not contain physical information. This is equivalent 
to saying that their values can be shifted arbitrarily by adding an appropriate 
counter-term. In our case this would correspond to shifting k. However, since this 
is a Chern-Simons term we only have the freedom to shift it by an integer amount. 
The fractional part of this term (in units of k) is therefore a physical observable 
[38]. It would be interesting to understand if this observable is constrained by the 
high-spin symmetry, like the other correlators discussed in [24]. 

Removing both contact terms, we are left with the parity-even result 

(J-(- 1 )J + ) = - N ^^. («) 
lb 7rA 

Note that this also changes sign at A = 1, consistent with our theory stopping to 
make sense (at least as a unitary theory) there. 

4.4 (jWjWjW) 

In this section we compute the 3-point function (j(°>(—q)J + (—q')J-(—q")), with 
q, q r , q" all in the 3-direction. We first need to compute the J+ vertex, 

(M-q)$(k)<l> i ) = S}V + (q,k ] \), (47) 

shown in figure 10. It can be evaluated using (30) and (43), and the result is 

V + (q, fc; A) = -L [lk 2 s + ql - q l e -^M2^ . (48) 

Returning to the 3-point function, the diagrams which contribute to it are 
shown in figure 11. In the diagrams that contain a seagull vertex, the final ra- 
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(i) (J) (k) 



Figure 11: Diagrams contributing to (J ( -°\—q)J+(—q')J-(—q"))- 
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dial integrals cannot be performed analytically. However, using the step-function 
relation 



Q(x - y)Q(y - z) + Q(y - x)G[x - z) = G(x - z)Q(y - z) , (49) 

one can see that these diagrams cancel in triples, 

(c) + (d) + (e) = (f) + (g) + (h) = (i) + (j) + (k) = . (50) 

The remaining diagrams (a), (b) are given by 

d 3 p V {q,p-q)V + {q , ,p)V 1 {q",p + q') 



2N 



N 



p 2 (p - q) 2 (p + q') 2 



vr 1 + e 



— iriX 



dp s e 



-2iX arctan(2j/)-2iA" arctan(2j/") 



H) Ylpl + q 2 + q' 2 + qq' 
{Ap 2 s + q 2 ){Ap 2 s + q" 2 ) 



+ 



N 



iq 



I _ e -7ri(A+A'+A") 



7T 1 + e - n 



iX 



AXq(q + q') 



(b) 



d 3 p d 3 k V {q,p)V 1 {q",k-q") (p + k) 3 



N 



n niX" 



7T 1 + e 



(2vr) 3 (2vr) 3 p 2 (p + q) 2 k 2 (k - q") 2 (p - k) + 

A e -2 *^ arctan(2y)-2iA" arctan(2j/") 

dp s 



q-q" 



-niX 



q" 



4p 2 s + q z 



Nq-q" 
7r AXqq" 



tan 



ttA 



(51) 



where now A' = Asign(g'),x' = njfny' = ^r, and similarly for q". In writing this 
we took A — > oo where possible to simplify the expressions. Since there are no 
divergences (as we shall see), this does not change the final result. 

The remaining radial integrals cannot be computed separately, but the integral 
in the sum (a) + (b) can be computed analytically. The final result is 



(J {0 \- q )J + (-q)J-(-q")) 



N 



,7TiA" 



iq 



1 _ g-7ri(A+A'+A") 



7T 1 



+ 



N 



-niX 



AXq(q + q') 
t(q + 2q') 



Nq-q" 
7T AXqq" 



tan 



ttA 

Y 



-TTi(X+X") 



(52) 



7T 1 + e -7r ^ 4:Xq'(q + q') 

By conformal invariance, this correlator should have one parity-even and one 
parity-odd structure [35]. Let us consider the A-even and A-odd parts separately. 



26 



4.4.1 Even Structure 



Let x = 2 ■ ^ ne ^" even P ar t of (52) is 



N 



' I A-even 



mmq)y) l- v 



8ir\qq'(q + q') cos(sign(g)x) 
- q' 2 sin(sign(g)x + 2sign(g')x) - q" 2 sin(sign(g)x + 2sign(g")x) } • 



By simple trigonometry we can write this as 



iVsin(TrA) 
8~ ttA 



q'q" 



+ 



qq" 



+ 



m 

qq 1 



+ 



sin 2 (^) tan(^) 1 



ttA 



m 



(53) 



(54) 



By conformal invariance we expect only a single parity-even structure, so we expect 
that the second piece is a contact term. Indeed, it is easy to check that 



w 



(55) 



is a conformally-invariant contact term (independent of the second momentum). 



4.4.2 Odd Structure 



Using similar methods, the odd part of (52) can be written as 



A-odd 



iN sin 2 (^) 



TTA 



1 1 

q' \q\ 



q' 



(56) 



In this case it is not important for us whether there is a contact term, because the 
A dependence of the odd structure (which is what we are interested in) cannot be 
changed by its presence. Indeed, this A dependence is determined by the last term 
in the square brackets, which cannot be a contact term. 



4.5 (J( 2 )J( 2 )) 

In this section we compute the correlator (T. (— q)T ++ ) with q = 0. We again 

introduce the notation of a squared vertex f/ 2 for T ++ , shown in figure 12 (similar 
to figure 9). The filled-circle vertex is the exact vertex shown in figure 

13. 
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U 2 {q,k;X) = % ® - J = X + 



k T. 



++ 



nr> + — — + 



+ d — x — 3 + 



5 reflections . 



+ — ^ 3 — + 



Figure 12: Diagrams of (T ++ (ft(f)) in which a gauge field is connected to T ++ . 
3 1 I 



Figure 13: Shorthand notation for the exact vertex (cjy <fiA + ) . 
We find that 

U 2 (q, k- A) = [3 (1 + A 2 ) k 2 + 3A%(A;3 + q 3 ) + 2i\ (2 + A 2 ) k s q 3 ] 5{ . (57) 

The 2-point function of the stress-energy tensor is then given by figure 8, replacing 

J_ with T and J + with T ++ . The exact T vertex is given by V 2 . Using (42) 

and (57), the result is 

(T__(-g)T ++ ) = JJL [3iqi (l - e 2darctan ( 2A '>) - 12Ag 2 A 

- 24i\ 2 q 3 A 2 + 16A(1 + 2A 2 )A 3 j . (58) 

The result should contain a single parity-even conformal structure, up to diver- 
gences and contact terms. The divergent pieces in (58) can be subtracted by 
counterterms of the form Ad IM g up d p g lJp , A 2 e pup g pa d p 'g u a , and A 3 g fiu g pu , where g^ v 
is the linearized background metric that couples to the stress-tensor. The finite 
A-odd piece, proportional to gf, can come from the contact term [38] 

(T pv {-q)T pa ) ~ [e ppX q\q u q a - q 2 b va ) + (ji u)] + (p «->• a) . (59) 
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The remaining finite A-even piece is 

iV|g 3 | 3 sin(7rA) 
128 ttA 

4.6 (j^j^j^y 



(60) 



In this section we compute the correlator (T (—q)J+(—q')J+(—q")), with all ex- 
ternal momenta in the 3-direction. We expect two parity-even conformal struc- 
tures, corresponding to the free bosonic and fermionic theories, and one parity-odd 
structure [35]. 

The diagrams are shown in figure 14. As in the case of (J^J^'J^), all 
diagrams that contain the seagull vertex cancel in triples. We are left with (a) + 
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(b) + (t/). After integrating over the p% and angular directions, they are given by 



N 

16V 



,2iA arctan 2A 



dp* 



(24p 2 + q 2 + q' 2 + q" 2 ' 



o 



Ap 2 s + q' 2 [l-e 



(b) 
00 



iN 

87r' 

iN 

87T* 



,2iA arctan 2k 1 ( Ji 



-2i\' arctan(2j/') 
-A 



(Ap 2 + q 2 )(Ap 2 s + q' 2 )(Ap 2 s + q" 2 ) 

2iX" arctan(2j/") 



4p 2 + g "2 (1 _ e - 



-2i\ arctan(2j/) 



(9" - 9) 



dp s 



2iX arctan 2h! ( J 



(q - q) 



4p2 + q 2 

A d Ps 
Ap 2 + q 2 



4\p s + iq' (l _ e -2iA'arctan( 2y ') 



3 — 2iA arctan(2?/) 



-2iA arctan(2j/) 



(61) 



The integral obtained by summing over (a), (b) and (b') can be computed 
analytically, and it has a linear divergence: Subtracting the divergence 15 we 
obtain the 3-point function, 

(T^(-q)J + (-q')J + (-q")) - ^ ' 



32tt A qq'q" 



- mX (q' 4 - q' 2 q" 2 ) + e~ mX ' (<f 4 - q' 2 q" 2 ) + 



q ,2 q „2 e - m (x>+x") _ e niX g 4 + q ' q " ^ + 7qfq >> + 4g "2) _ 12 \ 2 4 q " q 2 



(62) 



The A-even part can be written as 



N sin(TrA) 



A— even 



32vr A 

3 

+ 



VI3 



q'q" qq" 



-,"13 



qq' 



cos 



\q\ 



q'q" qq" 



/ 1 3 I „// 1 3 

+ ^4- 



-2 



q"\q'\ + q'\q"\ 



sm 



VTA 

2 

7rA 

T 



+ 



(63) 

qq' q J \ 2 / 

Here, the bosonic structure is the one multiplying cos 2 \Jy), since it is the one 
that survives when taking A — > 0. The fermionic structure multiplies sin 2 ( z y), up 
to possible contact terms. 
The A-odd part of (62) is 
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15 This divergence is subtracted by the counter-term g A + A + , which is related by £0(3) 
invariance to the term g^ A + A~ that we already used to subtract the divergence in (J + J_) 
(see equation (44)). However, our regularization preserves only SO(2) invariance in the x 1 — x 2 
plane (as well as parity duality) under which these terms are not related. The two subtractions 
are therefore independent. 
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The first term inside the brackets is a contact term, while the second term is the 
expected parity-odd structure. 



4.7 Correlators at the Critical Fixed Point 

In this section we consider the planar correlation functions in the "critical fixed 
point" of the bosonic vector model with Chern-Simons interactions. This fixed 
point is reached by starting with the theory we discussed above, turning on the 
relevant "double-trace" deformation ^(0V) 2 > an d flowing to the IR while tuning 
the IR scalar mass to zero [39, 40]. This is equivalent to adding an auxiliary field 
a with a a($(j)) coupling, or performing a Legendre transform with respect to the 
operator Alternatively, one can start with the usual "critical U(N) model" 
and couple it to the U(N)k Chern-Simons theory. 

In the planar limit the effects of the "double-trace" deformation (J^) 2 are 
rather simple. First, the scalar propagator receives corrections involving a chain 
of scalar loops, connected by the A 4 vertex, and ending in a tadpole (see figure 15). 
In the planar limit one can add gluon lines inside each scalar loop in this chain, 
such that the scalar loop has the topology of a disk in double-line notation. These 
corrections are all power-law divergent and independent of the scalar momentum. 
They are subtracted with a mass counter-term 0^0. 

Second, in correlators of gauge-invariant operators each insertion can be con- 
nected to a similar chain of scalar loops that ends on the rest of the diagram, as 
shown in figure 16. In momentum space the scalar loops factorize, and we can use 
our previous results to sum over the chains of scalar loops. For instance, using 
diagrammatics as above, or using conformal perturbation theory in the large N 
limit, the 2-point functions in the presence of the double-trace deformation can be 
written as 



conformal invariance. For s = s' = we can use the exact 2-point function of 




(65) 



The second line vanishes unless s — s' — 0, since ( J < -°- ) )a 4 =o = if s ^ from 



31 



Figure 15: Corrections to the scalar propagator in the critical theory, in the planar 
limit. Gluon lines can run inside the scalar loops. 




(34) to sum up the series, 
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To reach the IR fixed point we take A4 — > 00, so we expand in y- <C 1 and 
pick up the leading term (dropping the contact term N/X4). Defining the scalar 
operator of the critical fixed point to be = A 4 J < - ^ ) we obtain that 
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which is the 2-point function of a primary operator of dimension 2 as expected. 

Similarly, the 3-point functions with one scalar operator in the presence of the 
A4 deformation are given by 
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In particular using (54) (dropping the contact term) and (56) we obtain, 
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By similar methods, the correlator of three scalar operators at the critical fixed 
point evaluates to 

(j(°)(-g) J(°>(-g') j(°)(-g")) crit - oc (J^(-q) J^-^W)^! 
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The extra |g||g'||g"| factor cancels all the momentum dependence of (40) so we get 
a constant, which is a pure contact term. We therefore obtain that up to contact 
terms 



(J^\-q)J^\-q')J^(-q")) crit 
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For A = this is a standard result for the critical fixed point [41]. The vanishing for 
all values of A is consistent with similar results for the theory of fermions coupled 
to Chern-Simons gauge fields, and for the analogous computations in Vasiliev's 
theory of gravity [20]. 

All correlators that do not involve a insertion remain unchanged at the 
critical fixed point in the planar limit. 



5 Analysis of the Results 

Let us summarize the N and A dependence of the conformal structures that appear 
in the various correlators computed in section 4. 



2-point functions: 
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3-point functions: 
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The correlators (-)bos., (')fer. and (-) dd in the above equations were defined around 
equation (9) 16 . Note that we computed the correlators only for specific momenta, 



16 A correlator which involves J*- 1 ) in the free theory of a real field (boson or fermion) , is defined 
as 1/2 of the same correlator for a complex field. The vector current and energy-momentum 



34 



but using conformal invariance this is enough to determine the full correlators. 
One can verify that when A — > our correlators indeed go over to those of 2N free 
real bosons. 



5.1 Relation to the Results of Maldacena-Zhiboedov 

Let us map our microscopic couplings N, A to the parameters N, A of [24], by 
comparing our correlators to the ones we reviewed in section 2.1. In [24] the 
normalization of the energy-momentum tensor was chosen such that its 2-point 
function matches that of N free real scalar fields. Our stress-tensor is canonically 
normalized at any value of A, and since we can determine the overall normalization 
by requiring that in the A — > limit we have 2N free real scalar fields, we find 
from (75) 

N = 2N^^. (79) 

7rA 

The result (74) is then also consistent with (10), providing a consistency check on 
the assumption that the results of [ ] hold also for the odd-spin currents. Our 
other results will provide additional consistency checks for this assumption. We 
interpret N as the effective number of degrees of freedom in our theory, since (TT) 
is one way to define this number for a conformal theory; we find that it decreases 
as we increase the coupling, as expected, and that it goes to zero in the A — > 1 
limit (which is an infinite coupling limit using the Yang-Mills regularization to 
define the Cher n- Simons coupling). 

Next, we would like to compute A, by comparing our correlators to the ex- 
pressions of [24] for "quasi-boson" theories, written in (10), (11) and (13). Let us 
define the normalized correlator 

/ 7"(si) 7~( s 2) 7(S3)\ 

/j^j^j^a™ = w J _ J _ ; m 



Comparing our results to equations (10), (11) and (13) we have 
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tensor of a free complex fermion are denned as J M = ipj^ip, 7) t „ = |V'7( ( li ® v)ifi- With these 
normalizations the 2-point functions of J ^ and T^ v in the free boson and fermion theories agree. 
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Using (79), we then find that 
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The proportionality constant can be fixed, for instance, by requiring that A = 
corresponds to A = 0. The result is 
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where we have arbitrarily fixed the sign by assuming that A and A have the same 
sign. 

One can now compute several other normalized correlators to check the results 
(79) and (83). We find 
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All of these results are precisely consistent with the results of [24], using the values 
of N and A computed above. 

In the "quasi-boson" theory we also have the parameter 03 which is the co- 
efficient of the triple-trace term in d ■ This term affects only (J^J^J^). 
Expressing (34) and (40) in terms of A and N we obtain 
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From the 2-point function we see that to match with the normalization of [24] (see 
(11)) we must have = J^(i + A 2 ), and we then obtain from (86) 
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Therefore, matching to (13) we find 

N 2 V 87r 2 AV 16iV 2 cos 6 (f ) V Stt 2 / 
In particular, for A6 = 2Att 2 X 2 we get 03 = 0. 

We can similarly map our results of section 4.7 for the critical fixed point to the 
results of [24] for "quasi- fermion" theories (see the discussion below (14)). Since 
(TT) is the same in the critical and non-critical theories, N remains unchanged. 
To extract A q f of the quasi-fermionic theory consider the normalized correlator (see 
(67), (69)) 

Using (79) in the above equation gives the relation 

A * occos 2 ^). (90) 



Now, the A limit should correspond to A q f — > 00. This fixes the proportionality 
constant and we obtain 

A qf = cot f^j . (91) 

As we did above for the "quasi-boson" case, one can write down all the other 
normalized structures that we computed in section 4.7, and verify that they are 
all consistent with the results of [24]. 

So far we have matched the values of physical parameters to those of Maldacena 
and Zhiboedov using correlators at separated points. High-spin symmetry then 
determines all the 3-point functions at separated points. As we saw in section 
4.3, our theory also contains contact terms which contain physical information. 
It would be interesting to understand whether the high-spin symmetry constrains 
them, and to compare them to the fermionic theory. 



5.2 The Relation Between the Scalar and Fermionic The- 
ories 

It was shown in [24] that the correlation functions of the "quasi-boson" theory, 
which is equal to the free bosonic theory as A — > 0, become equal to those of 
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(the Legendre transform of) the free fermionic theory as A — > oo. Similarly, the 
correlators of the "quasi-fermion" theory (which is equal to the free fermion theory 
as A q f — » 0) become those of the critical O(N) (or U(N)) scalar model when 

A q f —> OO. 

Our results above imply that the limit of A — > oo corresponds to A — > 1, which 
is the maximal allowed coupling when regularizing the Chern-Simons theory with 
Yang-Mills terms. In this limit our results should thus correspond to a Legendre 
transform of the theory of iVf er . free complex fermions, for some value of iVf er .. Let 
us determine this value by matching (TT) between the two theories (note that in 
[24], the normalization is such that this correlator is the same for a free fermion 
and for a free boson). 

In the free fermion theory (and in its Legendre transform), we have simply 
(TT) = 2iVfer.(TT)i, where (TT)i is the result for a single real free boson or 
fermion. In our bosonic theory we need to take the limit A — > 1 while simultane- 
ously taking N — > oo. We can parameterize this limit by keeping &ym = k — N 
fixed while taking N — > oo. From (75) we then see that in this limit N — > 2&ym- 
Thus, we expect our theory to go over to the Legendre transform of the theory of 
-^Vfer. = &ym free complex fermions in this limit. It is nice to see that this result is 
always an integer, although since we derived it only in the large N limit with fixed 
N/k, it could be subject to shifts of order one (which must still give an integer). 

Next, let us determine the Chern-Simons level of this fermionic theory when 
we move slightly away from the free fermion point. In our bosonic theory, using 
(83) and taking the A — > 1 limit as defined above, we see that A ~ % (i-x) • ^ n 
the fermionic U(Nf eT .) theory at weak coupling, A q f was normalized in [24] so that 
A q f = f-jr 8 ^ (see equation (4.25) of [ ]). Since 3-point functions of operators 
with s > are independent of the Legendre transform in the fermionic theory, 
we can use this result also for the Legendre-transformed theory, but first we must 
translate from this "quasi-fermionic" variable A q f to the "quasi-bosonic" variable A 
used above. Matching the high-spin correlators (13), (14) (which are independent 
of the Legendre transform in the large N limit), we find that they are related 
by A = 1/A q f, and we then find that fcf er . = k (up to a possible sign 17 ). This 
implies that in the large N limit, the theory of N scalars coupled to a U(N)j~ 

17 This sign can be determined by comparing the signs of the one-loop corrections to the three- 
point functions of currents in the bosonic and fermionic theories, but we will not do this here. 
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Chern-Simons theory is equivalent to the (Legendre transform of the) theory of 
(k — N) fermions coupled to a - N)k Chern-Simons theory. If we translate 
the Chern-Simons level to the one defined using the Yang-Mills regularization, we 
obtain that the theory of iV scalars coupled to a U(N)k YM Chern-Simons theory is 
equivalent in the large N limit to the (Legendre transform of the) theory of kyM 
fermions coupled to a U(k YM ) N Chern-Simons theory. 

As a first consistency check on this statement, note that the Chern-Simons 
theories (without the matter) that we find on both sides of this relation are equiv- 
alent by level-rank duality [26, 27, 28] (this is true even at finite N, and certainly 
at large N). This is an important consistency check on the duality, since in the 
large N limit, the computations of many objects in these theories (like the S 3 - 
partition function or correlation functions of Wilson lines) are dominated by the 
Chern-Simons contributions, which scale as N 2 , and these must agree for the du- 
ality to make sense. Moreover, the level-rank duality of Chern-Simons theories 
exchanges Wilson lines in symmetric representations with those in anti-symmetric 
representations, which meshes well with the exchange of scalars with fermions. It 
would be interesting to see if one could perhaps derive the scalar-fermion duality 
by integrating out the scalars in one theory and the fermions in the other theory 
(at least at large N), expressing the results as correlation functions of Wilson lines 
in the pure Chern-Simons theory (see [42] and references therein), and seeing if 
these correlation functions are related by level-rank duality. It would be interest- 
ing to perform further tests of the duality, for instance by computing the effective 
potential on both sides, or by comparing the fractional contact term coefficients 
mentioned above. 

The duality described above may be viewed as a large N version of bosoniza- 
tion in three dimensions; a theory of fermions coupled to a Chern-Simons theory 
is described in a purely bosonic language (the theory of fermions without any cou- 
pling to Chern-Simons is described here as the infinite coupling limit of a bosonic 
theory). It is interesting to ask if this bosonization could be an exact equivalence 
also for finite values of N. Since the duality between the bosonic and fermionic 
theories exchanges weak and strong coupling, it is very difficult to test if this is 
true (given that we do not know how to perform exact computations at finite N) 18 . 

18 In the bosonic theory at finite N, or in the Legendre transform of the fermionic theory, the 
classically marginal coupling Ag has a non-trivial beta function and should be taken to its fixed 
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In two dimensional bosonization we know how to construct the fermion operators 
as solitons in the bosonic theory, and to prove the duality rigorously, but it is not 
clear how to do this in our case. In supersymmetric theories we can test similar 
dualities by comparing moduli spaces, chiral rings, and so on, but we do not have 
this privilege in our case. It is interesting to note that the transformation of the 
Chern-Simons group in our theory is the same as that of the Seiberg-like duality 
found in [30] for M = 2 supersymmetric Chern-Simons-matter theories (up to a 
shift by the number of flavors Nf which is not visible in our large N limit); of 
course in that case both sides contain both scalars and fermions coupled to the 
Chern-Simons theory, while in our case we have only scalars on one side and only 
fermions on the other side. It is interesting to ask 19 if the scalar- fermion duality 
could perhaps be derived by flowing from the supersymmetric duality, in which 
case we could confirm its validity for finite N (since the supersymmetric duality 
is believed to be valid also at finite N). 

One test that we can perform at finite N involves deformations of our theories. 
In this discussion, for simplicity, we use the definition of the Chern-Simons cou- 
plings using the Yang-Mills regularization. First, note that the fermionic theory 
we discussed cannot really be at level fcf er . = N, since in the presence of one flavor 
the level must be half-integer for the theory to be gauge-invariant [43, 44, 45]. 
Let us assume that the correct level is fcf er . = iV — |. Now, let us start with the 
fermionic theory coupled to a U(kYu)N-i/2 Chern-Simons theory, and deform it by 
a mass term to the fermions, M%p a %p a . At scales below M we can integrate out the 
fermions, and remain with a pure Chern-Simons theory, whose level depends on 
the sign of M [43, 44, 45]. For one sign we end up with a (topological) (7(^ym)jv 
Chern-Simons theory, and for the other sign with a U(kyu)N-i Chern-Simons 
theory. 

For the duality to be valid also at finite N, we must obtain the same low- 
energy theory also on the bosonic side. On that side we start with the critical 
bosonic theory, which can be viewed as the deformation of a theory of scalars 
coupled to U(N)k YM by (T(ft4>, and deform it by Ma. The auxiliary field a now 
serves as a Lagrange multiplier, enforcing 0^0 = —M, and we need to understand 

point. It was shown in [21] that such an IR-stable fixed point exists at large N and small A, and 
it would be interesting to understand exactly when it exists. 
19 We thank D. Kutasov for suggesting this. 
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the behavior of the bosonic theory with this constraint. For large values of iV 
and any value of A, one can show that for positive M the bosonic theory has a 
stable vacuum with unbroken U(N), in which the scalars are massive, so that at 
low energies we obtain the U(N)k YM pure Chern-Simons theory. For negative M 
there is no such vacuum, but there is an alternative vacuum in which one of the 
scalars condenses, and the gauge symmetry is broken to U(N — 1). In this case 
we find at weak coupling that the other (2iV — 1) gauge bosons become massive 
and dynamical (by swallowing scalar fields), as does the remaining real scalar, 
so we obtain at low energies the U(N — 1)& YM P ure Chern-Simons theory. These 
computations are done assuming that N is large, and in some cases also that the 
coupling constant is small, but since the low-energy theory is a topological theory 
labeled by discrete parameters, we expect to find the same low-energy theory for 
any iV and k. The two low-energy theories that we find here, for the two signs 
of M, are precisely equivalent (by level-rank duality) to the two theories that we 
found in the fermionic case, thus providing weak evidence for the validity of the 
bosonization also at finite N. 

It is natural to generalize to the case where we have Nf flavors of massless 
scalars/fermions in the fundamental representation of some U(N) group. The 
computations of the large N correlation functions that we computed above for 
this CcLSG £1X6 cL straightforward generalization of our computations in the previous 
sections, though now we have Nj operators of each spin, so it is not a priori 
obvious if the results of [2 I] can be applied 20 . In any case, the natural conjecture 
following from the discussion of the previous paragraph is that (using the Yang- 
Mills definition for the Chern-Simons coupling) the theory of Nf ■ N scalars coupled 
to the U(N)k YM Chern-Simons theory is equivalent (up to a Legendre transform) 
to the theory of Nf ■ &ym fermions coupled to the U(kyu)N-Nf/2 Chern-Simons 
theory. Translating back to our definition of the coupling, we map the theory of 
Nf scalars in the fundamental representation coupled to U(N)k to the theory of 
Nf fermions in the fundamental representation coupled to U(k — N)k-N*/2- This 
is very similar to the supersymmetric dualities of [30] . One difference is an overall 
shift in k by Nf/2, which is probably related to the one-loop contributions of the 

20 There are also new correlation functions that can appear when Nf > 1, such as terms 
proportional to f ABC in ( J^' A jW s J^ G ) (see [35]), and it would be interesting to compute 
them and to check if they are consistent with the duality. 
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fermionic fields to k. Another difference is taking Nf/2 — > Nf, which is related 
to the fact that in [30] Nf was the number of fundamental chiral multiplets, and 
also the number of anti-fundamental chiral multiplets, so the overall number of 
flavor fermions was doubled. A final difference is that in the supersymmetric case 
the sign of k changes under the duality; above we did not fix this sign, and it 
is plausible that also in the scalar-fermion duality one of the sides should have a 
negative Chern-Simons coupling (or equivalently, that if we keep positive Chern- 
Simons couplings, then the two sides are related by a parity transformation). 

Note that the results we present here are much stronger than most previous 
results on bosonization in three dimensions (see, for instance, [46, 47, 48, 49, 
50, 51, 52, 53]), which claimed that the low-energy limit of the theory of massive 
fermions coupling to a gauge field is given by a Chern-Simons theory (and had non- 
local bosonic actions at higher energies). In our case we claim that for massless 
scalars/fermions we have an exact equivalence of conformal field theories. For the 
special case of N = 1 and k = 1, our duality seems very similar to the duality 
studied in [54, 55, 56]; it would be interesting to understand this better, and to see 
if the methods of these papers can be used to study our duality more generally. 

5.3 Comments on the Thermal Free Energy 

As described above, the duality is completely consistent with all correlation func- 
tion computations done to date; preliminary computations of exact planar 2-point 
and 3-point correlators in the fermionic theory are also consistent with the duality 
presented in the previous subsection [ ]. However, the duality is not consistent 
(already at large N) with the form of the thermal free energy of the fermionic 
theory, computed in [20], as this form suggests a different relation between the 
bosonic and fermionic theories (as noted in [20, 24]). It is also not consistent with 
the thermal free energy of the scalar theory [57, 58] if one computes it by similar 
methods, as this gives a result that does not even go to zero when A — > 1. This is 
problematic independently of the duality, since we saw above that (T^ u T pa ) van- 
ishes in this limit, so the number of degrees of freedom should go to zero, and this 
is not visible in the naive computation of the thermal free energy. 

Our computations of the correlation functions pass many consistency checks; in 
particular they agree with the results of [24] which were computed by completely 
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different methods. On the other hand, there are so far no consistency checks for the 
thermal free energy computations. Thus, we claim that the existing computations 
of the thermal free energy are not correct. One possible problem is the light-cone 
gauge which these computations use. This is defined by an analytic continuation 
from Minkowski space, and it is plausible that such a continuation gives correct 
results for correlation functions (which are analytic in the momenta), but not for 
the Euclidean partition function compactified on a circle. 

Another possible problem involves the large N limit. The computations in 
question are performed in the limit of large N and large volume, and they take 
the fermions to be anti-periodic on the Euclidean thermal circle, and the scalars 
to be periodic on the circle, namely they assume that the holonomy of the Chern- 
Simons gauge field around the Euclidean thermal circle is trivial (Aq = 0). This 
assumption is expected to be correct at very high temperatures. It is also valid 
if we take the large volume limit first, so that the dynamics of the zero mode of 
the holonomy is decoupled. However, it is not clear if it is valid when we take the 
large N limit first. In particular, as discussed in [59], in vector models coupled 
to Chern-Simons theories, the holonomy becomes trivial (at least on S 2 ) only for 
very large temperatures T, obeying VT 2 ^> N (where V is the spatial volume). 
In the free theory this was explicitly checked in [59], where it was found that for 
temperatures that do not scale with N the holonomy is actually uniformly spread 
out on the thermal circle in the large N limit, leading to a vanishing free energy 
at order N. We expect this to remain true also for finite A. 

The computations of the thermal free energy in the Chern- Simons-matter the- 
ories in [20, 57, 58] use the standard 't Hooft limit, where only planar (disk) 
diagrams are kept. This is valid when we take the large N limit first, keeping 
everything else (like the volume or temperature) fixed, and take any other limits 
(like large volume) later. As described above, if we do this for the theory on S 2 
we land in the low-temperature phase, where the free energy vanishes at order 
N, so it does not give any useful comparisons between the bosonic and fermionic 
theories. (On other manifolds, like higher genus Riemann surfaces, the situation 
is more subtle since the Chern-Simons theory has many degenerate ground states 
[60], and we will not discuss this case here.) In order to get a non- vanishing free 
energy at order iV we need to be in the high temperature phase, with VT 2 ^> N, 
but then it is not obvious that the standard 't Hooft large N expansion applies. 
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This expansion is particularly subtle in our case, in which the leading order term, 
of order N 2 , is given by a topological theory. For instance, one may worry that for 
any A > connected diagrams with n scalar loops could scale as n powers of the 
volume, since in the pure Chern-Simons theory the correlators of n Wilson lines 
are independent of their positions, and this power of (VT 2 ) n could overcome the 
suppression by N 2 /N n in the high temperature phase 21 . It would be interesting 
to try to fix the thermal free energy computation, and to use it to test the duality. 
Similarly, it would be interesting if (as in supersymmetric theories [61]) one could 
find a way to compute exactly the free energy of our theories on a Euclidean S 3 (at 
least in the large N limit), as this could provide another useful test of the duality; 
this free energy for our theories at weak coupling may be found in [62]. 
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